Fracture Processes Observed with A Cryogenic Detector 
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In the early stages of running of the CRESST dark matter search using sapphire detectors at 
very low temperature, an unexpectedly high rate of signal pulses appeared. Their origin was finally 
traced to fracture events in the sapphire due to the very tight clamping of the detectors. During 
extensive runs the energy and time of each event was recorded, providing large data sets for such 
phenomena. We believe this is the first time the energy release in fracture has been directly and 
accurately measured on a microscopic event-by-event basis. The energy threshold corresponds to 
the breaking of only a few hundred covalent bonds, a sensitivity some orders of magnitude greater 
than that of previous technique. 

We report some features of the data, including energy distributions, waiting time distributions, 
autocorrelations and the Hurst exponent. The energy distribution appear to follow a power law, 
dN/dE oc E~ 13 , similar to the power law for earthquake magnitudes, and after appropriate transla- 
tion, with a similar exponent. In the time domain, the waiting time w or gap distribution between 
events has a power law behavior at small w and an exponential fall-off at large w, and can be fit 
oc w~ a e~ w ^ w ° . The autocorrelation function shows time correlations lasting for substantial parts of 
an hour. An asymmetry is found around large events, with higher count rates after, as opposed to 
before, the large event . 



INTRODUCTION 

In the Spring of 1999 preliminary runs of the CRESST 
dark matter search Q were carried out at the Gran Sasso 
Laboratory (LNGS), a deep underground laboratory for 
low background physics located in the Apennines. In 
these first runs of CRESST a phenomenon was observed 
which we believe may be of interest for the study of crack 
and fracture formation in brittle materials. CRESST is 
a cryogenic detector, working in the vicinity of 10 milli- 
Kelvin t 2j. In addition to being deep underground for 
shielding against cosmic rays, it is carefully designed to 
minimize effects of radioactive background. The detec- 
tor elements were large (262 gram) high quality single 
crystals of sapphire, with a strip of superconductor (W) 
evaporated on one surface to serve as a sensitive ther- 
mometer. This system, as shown by tests with gamma 
ray sources, detects single events in the sapphire with en- 
ergies in the range from about 1 keV to several hundred 
keV with good energy resolution ( 0.5 keV) and good 
time resolution (40 or 100 [is for the onset of a pulse). 

In order to reach these low temperatures it is impor- 
tant to eliminate the effects of any vibrations ("micro- 



phonics" ) that might deliver energy to the crystal. Thus 
in addition to special suspensions to isolate the appara- 
tus, the crystals are held very tightly in their holders to 
prevent any even microscopic frictional effects. In the 
data to be discussed here this was effected by small sap- 
phire balls held against the sapphire crystal by a plastic 
clamp. The plastic of the clamp, delrin, is known to con- 
tract substantially at low temperature, thus providing 
additional "tight holding" . An unanticipated result of 
the small contact area of the hard sapphire balls and the 
great force of the clamp turned out to be a cracking or 
fracturing of the sapphire. This was observed as follows. 

When the system was first brought into operation, an 
unexpectedly high rate of signal pulses was observed. Ini- 
tial fears that this might be due to an unexpected ra- 
dioactive contamination were relieved by the observation 
that even an unknown radioactive contamination must be 
Poisson distributed in time, while the unexpected pulses 
appeared rather to come in "bursts" or "avalanches" . 
Examination of the time distributions showed that they 
were indeed non-Poissonian. 

Pulse formation and fractures: The pulses themselves 
resembled those seen from good particle events. However, 
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this is a rather unspecific criterion, due to the operat- 
ing characteristics of the detector. There are essentially 
three steps in the production of a signal pulse 1) A rela- 
tively localized energy release within a short time, 2) A 
rapid degradation of this energy into a uniform "hot" 
(~ 10 °K) gas of phonons produced through phonon- 
phonon interaction and decay, as well as interaction with 
the crystal surface, 3) Absorption of the phonons in the 
thermometer strip. This leads to a heating with an 
increase of electrical resistance for the superconductor, 
which is finally read out by SQUID electronics. The re- 
sulting pulse shape is well described by a model employ- 
ing the various thermal and electrical parameters of the 
system As may be seen from this brief description, 
the pulse shape is essentially determined by the thermal 
responses of the system and not by the initiating event, 
as long as it is "fast" . Hence any release of a given en- 
ergy in the crystal in a short time (/x seconds) leads to the 
same pulse shape and so examination of the pulses does 
not lead to an identification of their origin. An extensive 
search for the origin of the pulses was finally successful 
when it was noticed that there appeared to be markings 
or scratches on the crystal at the contact points with the 
sapphire balls. When the sapphire balls were replaced by 
plastic stubs, which are evidently much softer, the event 
rate immediately dropped from some thousands per hour 
to the expected few per hour. 



These observations strongly suggest that the pulses 
were due to some kind of cracking or micro-fracturing 
phenomena in the sapphire crystal and/or its support 
balls. Indeed, examination under a microscope revealed 
a small crater with radiating irregular fissures extending 
sideways and down into the crystal. Damage to the sap- 
phire balls was also observed. Since the reduction in rate 
after the exchange of the sapphire balls was so large, we 
believe the data taken with the sapphire balls are essen- 
tially all fracture events. If we accept this crack or frac- 
ture hypothesis, our data then represent a large sample 
of well measured fracture events, under low background 
conditions, and with good time and energy determina- 
tion. 



Calibration runs: In order to calibrate the energy scale 
regular calibration runs were carried out. In these runs 
the system is left undisturbed and a radioactive source 
supplying 120 keV photons (which can penetrate to the 
detectors) is inserted in an external plug in the shielding. 
These photon-induced events can be selected by using the 
resulting 120 keV peak in the data. Since a radioactive 
source produces statistically independent events, that is 
Poisson statistics, these events provide a useful compar- 
ison when studying statistical properties of the data. 



ENERGY DISTRIBUTIONS 

We believe this is the first time that the energy re- 
lease in microfracture has been accurately measured on 
a microscopic event-by-event basis. 

It is to be emphasized that the cryogenic method pro- 
vides an absolute measurement of the total energy release 
in the fracture. This is to be contrasted with the study of 
acoustic emission in materials or seismic measurements of 
earthquakes. There the energy determination is necessar- 
ily indirect since there are various assumptions and un- 
certainties concerning production, propagation, and de- 
tection involved in translating the observed signals into 
the true energy of the event. On the other hand the cryo- 
genic method, essentially calorimetric in character, is a 
direct measurement of the full energy. The energy scale 
is fixed by the calibration with known sources and the 
resulting accuracy of the CRESST energy determination 
is on the order of a few percent yij. 

In addition to the directness of the energy measure- 
ment an important feature of the cryogenic method is 
its great sensitivity. The closest previous technique ap- 
pears to be the study of acoustic emission in materials 
There the smallest emitting region considered is 
on the order of a square micron Q . This will correspond 
to the breaking of ~ 10 7 bonds in the crystal. On the 
other hand, our energy threshold is typically some keV 
(Fig. 1). This corresponds to the breaking of only a few 
hundred or thousand bonds. Thus the cryogenic method 
appears to be many orders of magnitude more sensitive 
than previous technique. Small cryogenic devices can 
even be sensitive to energies in the eV range Q and it is 
possible that studies of this type involving stress release 
of only a few atoms are feasible Q • 

In Fig. 1 we show the differential distribution dN/dE 
for the number of events N per unit energy, for four data 
sets with two detectors from Run9. The straight line is 
the result of a power law fit 

dN/dE oc E-P (1) 

to the lowest curve, which yields (3 « 1.9. Similar results 
are found from fits to other data sets. From a total of 
seven sets examined (from Runs 9, 10 and 11) (3 ranged 
between 1.7 and 2.0. An interesting point is that the 
rates do not appear to differ greatly from one data set 
to another, despite the fact that different crystals and 
mountings are often involved. At 21 keV for example, 
the rates over the various data sets vary between 4 and 
11 pulses/keV-hr. 

A power law of this type, called the Gutenberg- Richter 
law |8|, is well known for the "magnitudes" of earth- 
quakes. Unfortunately the "magnitude" is a seismic am- 
plitude and not a direct measurement of the energy of an 
earthquake. Thus a simple comparison is not possible. 
However if one takes the prescription that the seismic 
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FIG. 1: Energy spectra from four data sets of Run9, with 53 
hr for the upper pair of curves and 28hr for the lower. The 
straight line shows a fit to the lowest curve oc E~P , yielding 
f3 « 1.9. 
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FIG. 2: Waiting time distributions. Upper curve: fractures, 
fit to oc w~ a e~ w / wa . Lower curve: photon-induced events 
from a calibration run, fit to oc e~ w ^ w ° . 



amplitude to approximately the 3/2 power |8(,|9| repre- 
sents the energy, and uses the power » 1.0 found for 
the integral distribution of earthquake magnitudes 0, 
it corresponds to/3~l + l~1.7, not far from our 
(3 w 1.7 — 2.0. Of course, the six orders of magnitude 
range available for seismic data is much greater than the 
one or two orders of magnitude available here. 

It should also be noted that such power law, that is 
scale free, distributions appear in many phenomena, of- 
ten related to an underlying fractal process [lOj. In the 
acoustic emission recordings of microfracture events in 
brittle materials for example, such a distribution typ- 
ically appears, with a somewhat lower exponent, ~ 



TIME SERIES 

Waiting Time Distributions: A useful quantity in the 
study of intermittent data such as the present is the 
"waiting time" w. To each event i we assign Wi, the 
time interval till the next event, and study the distri- 
bution of these intervals. Fig. 2 shows the waiting time 
distribution for detector 2 in a 28 hr data set of Run 
9. The distribution has power law behavior at small w 
and an exponential fall off at large w, and an accurate fit 
is obtained with dN/dw cx w~ a e~ w/w °, with a = 0.33. 
Similar results are found for other data sets with a in the 
range 0.25-0.5. The parameter wo determines the loca- 
tion of the crossover from power law to exponential and 
is essentially the inverse rate or average waiting time, 
with w = (1 — a)u'Q. Qualitatively similar results, with 
a near to or somewhat less than one, have been reported 
for earthquakes in California |ll) . 

For the simple case of Poisson statistics, one expects 



a waiting time distribution oc e ~ w / w o ^ where 1/wq is the 
average count rate. The lower curve of Fig. 2 shows the 
waiting time distribution for the photon-induced events 
of a calibration run, with a fit to oc e ~ w / w o _ As expected 
there is a good fit, and with 1/wq in agreement with the 
event rate. 

An interesting point concerns the behavior of wo for 
fracture events as the the energy threshold for the sam- 
ple is raised. It appears that the form w~ a e~ w ^ w ° is 
preserved, with a varying little. Since the count rate is 
reduced however, the value of wq increases and so the 
crossover between power law and exponential behavior 
moves to larger w. Indeed, taking a given data set (Run9- 
d2,100/zs), repeatedly raising the energy threshold and 
fitting for w , we find a linear relation between the in- 
verse count rate, that is w, and the fitted wq. The slope 
and the relation w = (1 — a)wo then gives a global de- 
termination a « 0.26. 

The power law behavior for the waiting times at small 
w, as well as the power law for the energy distribution in 
the previous section, is suggestive of an underlying scale- 
free processes without any intrinsic dimensional parame- 
ter, as is common in fractal processes |l(i|. However this 
cannot be entirely true here since wq is a time and has di- 
mensions. Since e~ w / w ° corresponds in fact to a Poisson 
distribution, this may suggest an interpretation in terms 
of some basic scale free processes where several such pro- 
cesses are occurring independently and simultaneously 
and so are overlapping in the data. This arises trivially 
if the signals originate from more than one of the support 
points of the crystal, of which there were several; but one 
can also imagine independent crack systems beneath one 
support point. 

The increase of wq as the count rate goes down sug- 
gests that the limit of zero count rate is a kind of critical 
point: the waiting time becomes infinite as the distribu- 
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FIG. 3: Autocorrelation function C for the event rate from 
calibration data of RunlO. The lower curve is for photon- 
induced events (events in the 120 keV peak), and the upper 
curve is for fractures plus some admixture of compton scat- 
ters (events below the peak). For the photons the data is 
consistent with C = for (t — t') ^ as expected for Poisson 
statistics, with C(0) = Variance = R. 



tion becomes non-integrable and completely scale free, 
while 1/wo appears as a diverging correlation length. Un- 
derstanding wq is an interesting point for further study. 

Correlations in Time: We expect the existence of 
correlations in time, corresponding to the "bursts" or 
"avalanches". We use the event rate Rt of a calibration 
run to construct the autocorrelation function 



C(t-t') = (R t - R)(R? — R) 



(2) 



and compare C for photon-induced events and fractures 
in Fig. 3. While for photons we have C = as expected, 
for the microfractures there are correlations lasting for 
substantial fractions of an hour. These long-term cor- 
relations are found for the fracture events of all data 
sets. The physical origin of the correlations may be in 
stress relaxation phenomena where a slow "diffusion" of 
strain ^1 can trigger new microfractures when meeting 
other weak spots in the crystal. 

Hurst exponent: The autocorrelations as in Fig. 3 can 
be approximately fit to power laws oc (t — t')~ p . As noted 
above, this is suggestive of the scale free, self-similar be- 
havior associated with fractal statistics. A way of char- 
acterizing such behavior is in terms of what is called the 
Hurst exponent H; and we can check the plausibility of 
such a description by comparing the consistency of H 
found in different ways. Table I shows H found in three 
ways for various data sets. First the autocorrelation ex- 
ponent p is fitted to find H = 1 —p/2. The next column 
shows H determined by the "growth of the standard devi- 
ation" , a characterization of the fluctuations in the event 
rate ~ t , where t ' 5 would be the classical gaussian or 



TABLE I: Exponent H found by different methods, dl and 
d2 refer to the two detectors in operation, and 40, 100 /j,s to 
different digitization windows used in data taking in Run9. 



Data Set 


Autocorr. 


Stnd. Dvtn. 


Sh. Entropy 


Run9 dl W0[is 


0.77 


0.70 


0.69 


Run9 d2 100^s 


0.80 


0.80 


0.80 


Run9 dl 40^s 


0.73 


0.70 


0.67 


Run9 d2 40/is 


0.69 


0.70 


0.65 


RunlO d2 


0.59 


0.63 


0.59 


Runll dl 


0.60 


0.64 


0.53 


Run 11 d2 


0.69 


0.66 


0.62 



random walk behavior with finite range correlations. Fi- 
nally, the last column gives H found from the "Shannon 
Entropy" , related to the probability of the number of 
events over a time interval t [13j • Although the fits were 
not all excellent and there is considerable fluctuation in 
the results, the overall rough consistency of the three 
determinations supports the picture of a scale free, self- 
similar process. We do not necessarily expect the same H 
for different data sets since these involve different energy 
thresholds and sensitivities. 



CLUSTERS 

A frequently used concept in the earthquake literature 
is the "Omori Cluster": a "big shock" followed by "af- 
tershocks". As Fig. 1 shows, and as is also the case for 
earthquakes, there is no separate class of high energy 
events — no distinctive "big shocks" . Naturally, as should 
be expected from the "avalanches" or correlations, given 
any event, there is a general increase in rate at nearby 
times. Although this increase is quite substantial, (a fac- 
tor four with one second bins, see Fig. 4) this simply 
reflects the "bursts" or "avalanches" and is not specific 
to "big events" . 

More specific to "big events" , however, we find a time 
asymmetry with respect to "before" and "after". That 
is, there is on average more activity after, as opposed to 
before, "big events" . Fig. 4 shows the count rate from a 
data set of Run9, for times close to "big events" , with "af- 
ter" (upper histogram) and "before" (lower histogram) 
plotted separately. "Big" was defined as a pulse with 
E > 300fceU. The bin size is 2.0 x l(T 4 hr= 0.72s. One 
notes a significantly higher rate in the first bin "after" 
relative to that in the first bin "before" . There is a rea- 
sonable fit to a power law for the decline in the rate to- 
ward the average value, and a significantly steeper power 
"after" relative to that "before" is found. Similar results 
are obtained for other data sets. 

An asymmetry of this type appears to exist in seismic 
data and certain models |l4| and seems to indicate that 
the "big events" tend to occur early in the "bursts" . 
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FIG. 4: Count rates in the vicinity of "big events" , showing a 
time asymmetry "before" and "after" the "big events" . From 
a 53 hr data set of Run9, plotted in 0.72 s bins. The upper 
(dotted, red) histogram is for times following the "big event" 
and the lower (solid, blue) histogram for times preceding the 
"big event". The straight lines are power law fits, yielding a 
power 0.13 ± 0.01 "before" and a power 0.25 ± 0.01 "after". 
There were 1082 "big events", defined as a single pulse with 
E > 300fceV\ The average rate in the run was 526/hr. 



CRACK PROPAGATION AND MATERIAL 
PROPERTIES 



Our material is a single crystal of high purity [15(. In 
crack propagation models the growing stress enhance- 
ment at the crack tip implies that a "hard spot" is nec- 
essary to limit the propagation of a crack; thus when 
a homogeneous stress is applied to a defect-free material 
there is nothing to stop a propagating crack. Presumably 
the microfractures here were limited by the random, non- 
homogeneous stress and defect field which quickly arises 
as fractures form in the pure material. This may have 
been assisted by the damage to the small sapphire balls, 
leading to an irregular application of the stress. Although 
we speak of "cracks" , it should be kept in mind that from 
our simple observation of pulses we cannot infer the ex- 
act nature of the microfracture. Finally, with respect to 
materials it should be noted that our system is of course 
quite opposite to those in the geological context, where 
one has highly heterogeneous systems, while here we have 
a very pure material. 



DEVELOPMENT OF THE TECHNOLOGY 

It is interesting to contemplate the extension of this 
method in the study of fracture phenomena. The su- 
perconducting thermometer, and perhaps other cryo- 
sensors Q, can be applied to many materials. The very 
low temperature and large crystals of the dark matter 
search would not always be needed, and indeed it might 
be possible to follow the crack development in time with a 
smaller and thus faster system. However, low background 
conditions may still be necessary to avoid contamination 
of the data by non-fracture events. In the present data 
the crystal was contacted by several of the small sap- 
phire balls, and we are unable to determine where an 
event originates. Such effects lead to a dilution of cor- 
relations, which thus may be intrinsically much stronger 
than appear here. In an apparatus especially designed for 
such studies one could arrange to have only one "hard" 
contact and with a known force. Finally, since the energy 
range available is relatively small compared to that for 
earthquakes it would be useful to consider techniques for 
increasing the dynamic range. 
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